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Abstract Motivated by the wide applications of thermal stability and phase transition, we investigate 
thermodynamic properties of charged BTZ black holes. We apply the standard method to calculate the 
heat capacity and the Hessian matrix and find that thermal stability of charged BTZ solutions depends 
on the choice of ensemble. To overcome this problem, we take into account cosmological constant as a 
thermodynamical variable. By this modification, we show that the ensemble dependency is eliminated 
and thermal stability conditions are the same in both ensembles. Then, we generalize our solutions to 
the case of nonlinear electrodynamics. We show how nonlinear matter field modifies the geometrical 
behavior of the metric function. We also study phase transition and thermal stability of these black 
holes in context of both canonical and grand canonical ensembles. We show that by considering the 
cosmological constant as a thermodynamical variable and modifying the Hessian matrix, the ensemble 
dependency of thermal stability will be eliminated. 


1 Introduction 

Discovery of the three dimensional BTZ (Banados-Teitelboim-Zanelli) black holes is accounted as one 
of the greatest achievements in gravitational framework T|. These solutions provide simplified machin¬ 
ery that enables us to investigate different aspects of the black objects, such as their thermodynamics 
mm- The generalization of these black holes to higher dimensions has been employed to investigate 
the black holes solutions and their thermodynamics [5]. In addition, BTZ black holes have been used 
to expand our understanding of gravitational interaction in low dimensional spacetimes [B]. In context 
of the string theory and quantum gravity, there have been several studies regarding the effective action 
of the string theory and BTZ black holes [7]. Due to the fact that these solutions being asymptot¬ 
ically AdS, there has been various studies regarding AdS/CFT correspondence [5]. In addition, the 
noncommutative geometry in three dimensions has been of interest recently jB] and it was shown that 
gravitational Aharonov-Bohm effect could be originated from noncommutative BTZ black holes |10j . 
Moreover, in order to obtain quantum aspect of three dimensional gravity, entanglement and quantum 
loop entropy of the BTZ black holes have been investigated in literatures m- Also, the holography of 
the BTZ black holes has been investigated in details [T2I . 

Nonlinear theories have been of interest due to their special properties. Most of natural systems have 
nonlinear behavior and usually governed by nonlinear theories. Among these theories, the nonlinear 
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electrodynamics, specially Born-Infeld (BI) types, have been of interest [TJ:. BI electrodynamics first 
has been introduced in order to solve the self energy of charged point like particles in Maxwell theory 
M- It was also shown that this theory enjoys the absence of the shockwave m ■ On the other hand, it 
was also proven that there is no birefringence phenomena in this nonlinear theory of electrodynamics 
[161 . In addition, the electric-magnetic duality is another important property of this nonlinear theory 
m Another important property of the BI types of electromagnetic field is in context of the string 
theory. It was shown that the Lagrangian of these theories may be arisen in limit of low energy effective 
of heterotic string theory Q]JJ. One of the interesting properties of the BI types theories is the fact that 
their series expansion for large values of the nonlinearity parameter yields the same structure; The first 
two terms of their series expansions are Maxwell invariant and a quadratic Maxwell invariant coupled 
with first order nonlinearity parameter with different coefficients m■ Therefore, one can construct 
a generalized theory of nonlinear electromagnetic field consisting the quadratic Maxwell invariant in 
addition to the Maxwell Lagrangian [20]. Besides, it is arguable that in generalizing linear Maxwell 
theory to nonlinear one, we should take only small values of nonlinearity into consideration. Also, in 
context of experimental results, only small power of nonlinearity are applicable. Another interesting 
property of this theory is the fact that we are dealing with nonlinearity as a perturbation. In other 
words, we are considering additional term as a perturbation to the Maxwell Lagrangian. Therefore, we 
are dealing not only with a nonlinear theory but also a perturbative one too. 

The concept of black holes as thermodynamical systems enables one to consider various thermo¬ 
dynamical aspects of solutions. One of the most important thermodynamical properties of the black 
holes is their thermal stability. For black holes being physical objects, they should be stable in con¬ 
text of dynamical and thermodynamical frameworks. The instability of black holes means whether the 
system is completely non physical or it may have phase transition. In other words, the system may go 
under phase transition to acquire a stable state. The stability of BTZ black holes and black string have 
been studied in several articles [2Tj . Myung showed that there is a possible phase transition between 
non-rotating BTZ black holes and a thermal AdS spacetime [22]. In Ref. [22. Hawking-Page phase 
transition of BTZ black holes and thermal soliton were studied. Regarding the BTZ black holes with 
torsion, it was shown that phase transition depends on theory of gravity under consideration [2T. 

The rest of the paper is organized as follows. In the next section, we present a brief review of 
the Lagrangian and related field equations of Einstein-Maxwell gravity as well as Einstein-nonlinear 
Maxwell theory. In Sec. [3] we investigate thermodynamics properties of charged BTZ black holes and 
discuss ensemble dependency. Then, we present a suggestion for overcoming the mentioned dependency. 
We generalize our results in Sec.[4]to the case of nonlinear electromagnetic field. It means that we study 
thermal stability of these black holes in context of both canonical and grand canonical ensembles, and 
suggest that in order to avoid ensemble dependency, we should consider the cosmological constant as 
a thermodynamical variable. We finish our paper with the highlight of conclusions. 


2 General formalism and field equations 

In order to study 3-dimensional black holes in the presence of a matter field, we employ the following 
Lagrangian 

Ltot = R — 2 A + L m , (1) 

where R and A are, respectively, the Ricci scalar and the (negative) cosmological constant. The last 
term in Eq. m is the Lagrangian of matter field, which we choose an electromagnetic field. The usual 
linear electrodynamics is Maxwell field with the following Lagrangian 


L(J-) = -F, (2) 

where the Maxwell invariant is F = F a bF ab in which F a & = d a Ai, — di,A a is the electromagnetic 
field tensor and A& is the gauge potential. In addition to the linear Maxwell field, one can use the 
nonlinear models. Nonlinear models of electrodynamics were introduced with different motivations. In 
this paper, we consider a quadratic correction in addition to the Maxwell Lagrangian to obtain the 
following nonlinear electrodynamics 


L(F) = -F + pF 2 + 0(f3 2 ), 


( 3 ) 
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where j3 is nonlinearity parameter. As one can see for vanishing nonlinearity parameter, Maxwell 
theory of electromagnetic field is recovered. Using variational principle and varying Lagrangian © 
with respect to metric tensor and gauge potential, we can find the following field equations 

Rab ~^^9ab F Ag a b = m ^9abL / (F') ‘IL j^F ac F^ , (4) 

da {V~gL^F ab ) = o, (5) 

where Ljr is derivation with respect to Maxwell invariant. In order to obtain 3-dinrensional solutions, 
one can employ the following static metric ansatz 

ds 2 = -f( r )dt 2 + ^- + r 2 de 2 . (6) 

.fin 

Regarding the structure of the electromagnetic field, we are interested in radial electric field. There¬ 
fore, we consider the gauge potential as 


A a = h(r)SS, (7) 

which results into the following nonzero components of electromagnetic field 

F tr (r) = -F rt (r). (8) 


3 Thermodynamic properties of charged BTZ black hole 

As we mentioned in Sec. [Q static BTZ black hole solutions in the presence of Maxwell field have been 
obtained before 1] . In addition, thermodynamic properties [2] and heat capacity [3] of charged BTZ 
black holes have been investigated in literature. In other words, thermal stability of the solutions has 
been investigated in canonical ensemble. In this section, we review the BTZ black hole solutions in 
the presence of Maxwell theory. Then, we show that there is a case of ensemble dependency in study¬ 
ing thermal stability of the solutions which will be removed by special consideration of cosmological 
constant as a thermodynamical variable. 

In order to obtain electric potential, one should use Eqs. © and 0 with considering the linear La¬ 
grangian of Maxwell theory, © and 3-dimensional metric ®. These considerations yield the following 
forms for the electromagnetic potential and the electric field 


h{r) = gin (y) , 

(9) 

F tr = E{r) = —, 

(10) 


in which q is integration constant related to total charge of black holes. It is notable that we inserted 
a scale factor, l, into Eq. © to obtain dimensionless argument for the logarithmic function. Although 
this scale parameter does not affect the consistency of field equations, we will find that it has a decisive 
rule for removing ensemble dependency. Considering obtained electromagnetic tensor, one can show 
that tt and rr components of Eq. © yield the same differential equations. We find 


= jtt = r 2 f + 2 Ar 3 + 2 q 2 r = 0, 

(11) 

jee = r 2 f" + 2Ar 2 - 2 q 2 = 0, 

(12) 


where prime and double prime denote first and second derivative with respect to r, respectively. It is 
notable that jgg = j' tt — ^ j t t , and therefore, the solutions of Eq. (flTTl satisfy Eq. ED. simultaneously. 
The consistent solution is [TJ 

f(r) = —Ar 2 - 2 q 2 In (y) - to, 

where the integration constant m is the geometrical mass of black holes. 


( 13 ) 
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Fig. 1 f(r) versus r for A = —1, l = 1, q * 1, and m = 0.1 (solid line), m = 1 (dotted line) and m = 3 
(dashed line). 

In order to interpret the solution as a black hole, we should obtain its singularity(ies) and horizon(s). 
For investigating the existence of the singularity, we obtain the Kretschmann and Ricci scalars which 
yield the following relations 


= 12yl 2 + ^ + 


R = 6/1 + 


2 q 2 


(14) 

(15) 


Taking into account Eqs. m and m, it is evident that there is an essential singularity at the origin. 
Besides, considering the behavior of the Kretschmann and Ricci scalars at large values of r, one can 
confirm that obtained solution is asymptotically AdS. Now we should discuss the existence of horizon. 
Plotting metric function versus r, we find that the number of the metric function’s roots (horizons) is 
a function of free parameters (see Fig. IT}). Fixing l and q , we find that for sufficiently small geometrical 
mass, there is no root and obtained solutions are presenting naked singularity. There is an extremal 
mass m ex t in which only one extreme root is observed. Increasing this value lead to a black hole with 
two horizons (one inner (Cauchy) horizon and one outer (event) horizon). 

In order to find the values of the horizon radius, one should consider the vanishing metric function 
at the horizon. Therefore, for the case of / (r+) = 0, one can find the following relation 


m + Ar± + 2 q 2 In 



= 0 , 


(16) 


in which r_ and r + are the horizon radii with the following explicit forms 


r± = l exp 


TO L\y± 

' V 2~ 


(17) 


(-?) 


and L\y- = LambertW 


-1 


Ar 


exp 


(-?)] ( fo 


Dr more 


where Lw+ = LambertW exp | 

details about LambertW(x) function see Ref. |l25l). 

Now, we are in position to investigate the thermodynamical properties of the solution and study its 
thermal stability in context of both canonical ensemble (heat capacity) and grand canonical ensemble 
(Hessian matrix). 
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3.1 Thermodynamics and conserved quantities 

Using area law for calculating the entropy of the system in Einstein gravity US], leads to the following 
result 


S = 


irr + 

~ 2 ~ 


(18) 


On the other hand, considering the Gauss’s law for obtaining total electric charge leads to 


Q = I’ 


(19) 


while for the electric potential, one can use the following standard relation 


^ = A^ 


Preference 




( 20 ) 


In order to calculate the temperature of the black hole, we use the definition of the surface gravity, 
which in the case of our static black hole will be [2] 


/' (r+) = Ar 2 + + q 2 
4-7T 2nr + 


( 21 ) 


Due to our interest in AdS solutions and inspired by AdS/CFT correspondence, we use the coun¬ 
terterm method in order to calculate finite total mass of black hole. As one can see, evaluating action of 
this configuration at infinity, yields infinite value. In order to overcome this problem, we add additional 
terms to action. One can show that for the obtained solution with flat boundary, R a bcd{ 7 ) = 0, the 
finite action is 

dfinite = Ig + lb + let , (22) 

where the bulk, boundary and counterterm actions are, respectively, 

Ig = r- [ d 3 xy/=gL tot , (23) 

lOTT J M 

h = f d 2 x^K, (24) 

JdM 

!ot = d 2 xy/^y (y~A) , (25) 

which 7 and K are, respectively, the traces of the induced metric, j a bi and the extrinsic curvature, 
K ab , on the boundary dM. Using the Brown-York method of quasilocal definition with Eq. m-m, 
one can introduce divergence-free stress-energy tensor as follow 


T 


ab 



K~i ab + 


(26) 


Then, the quasilocal conserved quantities associated with the stress tensor of Eq. (ESI) can be written 
as 

Q(0 = j de^T ab n a i b : (27) 

B 

where n a is the timelike unit normal vector to the boundary B. Taking into account Eq. m with 
£ = d/dt as a Killing vector, we find following result for obtained solution 



( 28 ) 







6 


3.2 Thermal stability of charged BTZ black hole 


In order to investigate thermal stability of the black holes, one can adopt two different approaches to 
the matter at hand. In one, the electric charge is considered as a fixed parameter and heat capacity 
of the black hole will be calculated. The positivity of the heat capacity is sufficient to ensure the local 
thermal stability of the solutions. This approach is known as canonical ensemble. In this case, the 
system which is unstable may go under phase transition to stabilize. The phase transition points are 
where the heat capacity has root(s) or diverges. Another approach for studying thermal stability of 
the black holes is grand canonical ensemble. In this approach, the thermal stability is investigated 
by calculating the determinant of Hessian matrix of M(S,Q ) with respect to its extensive variables. 
The positivity of this determinant also represents the local stability of the solutions. In what follows, 
we study stability of the solution in context of both ensembles. Before we conduct our study, it is 
worthwhile to mention an important point. In some cases, these two approaches may admit the stable 
phase for black holes. But studying the positivity of these two quantities is not sufficient for ensuring 
the physical thermodynamical behavior of the system. In order to have more realistic results and 
also enriching them, studying the behavior of the temperature is necessary, simultaneously. In other 
words, positivity of the temperature may denote the systems being physical whereas its negativity is 
representing non-physical systems. 

For the canonical ensemble, the heat capacity is Cq = T (§§*)q, where by using chain rule, one 
can rewrite it in the following form 

(m.) 

\ dr +Jn 

C Q = T-. -(29) 

( MIL] 

\ dr +J Q 

It is a matter of calculation to show that the heat capacity for charged BTZ black hole will be J3] 


7rr_|_ (Ar\ + q 2 ) 
2 (A4 - q>) 


(30) 


It is evident from the obtained relation for heat capacity that, in case of the asymptotically AdS 
spacetime (A < 0), the denominator of Cq is negative. Therefore, in order to have positive heat 
capacity (stable charged BTZ black hole), the following inequality must be hold 

(J 2 

A<-\. (31) 

r + 


In order to have type one phase transition (Cq = 0), we obtain the following phase transition point 

n + - ^ (32) 

In addition, we should note that only in case of asymptotically AdS spacetime the phase transition 
type one exists. As for the phase transition type two (Cq —> oo), one can find divergence points of the 
heat capacity by the following relation 

r2+ = vX (33) 

To summarize, we interestingly find that the phase transitions for asymptotically AdS or dS solu¬ 
tions are, respectively, the first type (vanishing point of heat capacity) or second type (singular point 
of heat capacity). In other words, there is no divergency for the heat capacity in case of asymptotically 
AdS spacetime. 

Now, we plot various diagrams in order to investigate thermal stability of charged BTZ solution 
in context of heat capacity. Fig. [2] shows that there is a critical value (ri+) in which for r + < n+ 
the temperature (and the heat capacity) is negative, and therefore, the black hole is non-physical. 
Otherwise, the system has positive heat capacity and temperature. In other words, the black hole is 
physical and in thermally stable phase for r + > ri+. 
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Fig. 2 Cq and T (bold lines) versus r+ for A = — 1, and q = 1 (solid line), q = 1.2 (dotted line) and q = 1.4 
(dashed line). 


Now, we study the stability of solutions in the case of grand canonical ensemble. In this case, the 
standard extensive parameters are charge and entropy. Therefore, one can see the mass of the black 
hole will be in the following form 


M(S,Q) = -- 


AS 2 




and one can obtain related Hessian matrix with the following form 


IT M _ 

“S.Q — 


which by using Eqs. (US and ra its determinant will be 


r -AS 2 +Q 2 ir 2 

2 Q ~\ 


/ S \ 

2 Q 

S 

Ki¥)J 


H ' M I _ 

s,q\ — 


(Ar\ - q 2 ) . / r 


■n 2 r\ 


‘“'if'- 


4 q 2 
7r 2 r)_ 


(34) 


(35) 


(36) 


Now, we plot Fig. [3] for studying thermal stability in context of grand canonical ensemble (Eq. 

M- 

It is evident that there is a case of ensemble dependency. In other words, considering the same 
parameters, these two approaches for studying thermal stability do not yield the same result. 

The ensemble dependency may be originated from one of the following reasons: first of all, because 
of the dependency of the Hessian matrix to thermodynamical variables, our choices of thermodynamical 
variables may be wrong. In other words, total mass of the black hole may depend on more extensive 
parameters that we had taken into account. Second, the system may be constructed in a way that being 
ensemble dependent is a part of its fundamental properties. In this paper, we take the first assumption 
into account and try to solve the ensemble dependency problem. 

Motivated by recent analogy in which (negative) cosmological constant is considered as a thermo¬ 
dynamical variable, we take cosmological constant as an extensive parameter. Therefore, the Hessian 
matrix from 2x2 changes into a 3 x 3 one. As we pointed out, we inserted the scale factor, l, into 
the solutions for obtaining dimensionless argument for the logarithmic function and in general it is 
independent of cosmological constant. But since dimension of cosmological constant is ( length )~ 2 and 
in order to overcome ensemble dependency, we use the following logical relation for the scale factor 

r 1 = V^A. (37) 















Fig. 3 Different scales of |H^q| and T (bold lines) versus r+ for l = 1, A = — 1, and q = 1 (solid line), q = 1.2 
(dotted line) and q = 1.4 (dashed line). 


By this consideration, the structure of Hessian matrix for our study will be 


ttM 

n S,Q,A 


—AS 2 +Q 2 -n 2 

7T 2 S 2 

2 Q 
S 

s 


2 Q 
' S 


In ( — 


^ _ 


4 AS 2 ) 
A 


Q 
A 

_2l 

2 A 2 J 


with the following determinant 


IttM I 


[(2 Ar\ — q 2 ) In (—Hr 2 ) — 6g 2 ] (Hr^_ + q 2 ) 


87r 2 yl 2 r^_ 


(38) 


(39) 


Now, in order to study the behavior of Eq. (1551) . we plot Fig. 2) It is clearly seen that the ensemble 
dependency is solved by this consideration. It is true that one may argue that for small values of the 
horizon radius there is a region in which the determinant of the hessian matrix is positive, hence black 
hole is stable which is in contradiction to earlier result that was derived for heat capacity. It is notable 
that this region is located at the place where the system has negative temperature. Therefore, we take 
away this region for the reason of non-physical black hole solution. 


4 Generalization to nonlinear electrodynamics 

In this section, we consider the generalization of linear electrodynamics to nonlinear one by adding 
the quadratic Maxwell invariant to the Maxwell Lagrangian. Using Eqs. © and 0 with Lagrangian 
of nonlinear electromagnetic field 0 and metric 0, one finds 

h(r) = qln (y) + + O (/? 2 ) , (40) 

in which q is integration constant related to total electric charge of the solution. Using obtained function 
for the electric potential results into electric field F tr to be 

F tr = E(r) = q --^-l3 + 0(p), 


(41) 
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Fig. 4 Different scales of 10 2 1 H sf< 3 ,/i| and 10T (bold lines) versus r+ for A = —l 2 = — 1, and q = 1 (solid 
line), q = 1.2 (dotted line) and q = 1.4 (dashed line). 


which for small values of /3, obtained relations reduce to linear Maxwell case. 

Now, we are in a position to obtain metric function of this nonlinear theory up to first order of 
nonlinearity parameter. To do so, we use Eq. 0 by applying electromagnetic field tensor m which 
will lead to the following differential equations 

e rr = e tt = r 2 f + 2Ar 3 + 2 q 2 r - -y^/3 + O (/ 3 2 ) , (42) 

egg = r 2 /" + 2Ar 2 - 2 q 2 + ^/? + O (/3 2 ) . (43) 

It is easy to show that egg = e' tt — yejt, and therefore, it is sufficient to solve Eq. (THZl) . One can show 
that metric function for this configuration will be 

f(r) = -m - Ar 2 - 2 q 2 In (y) - ^-/3 + O (/ 3 2 ) , (44) 


where m is an integration constant which is related to total mass. 

The next step will be devoted to study the possibility of existence of singularity and horizon for 
the solution. In other words, obtained solution may be interpreted as a black hole if there exists a 
curvature singularity covered with horizon. In order to investigate the existence of singularity, one can 
study curvature scalars such as Kretschmann and Ricci scalars. Studying these scalars enable us to 
investigate the asymptotical behavior of the solution too. Evaluating Kretschmann and Ricci scalars 
result into 


« = 6 ' 1 + |r + |r' s + 0 (/ j2 ). 




(45) 

(46) 


which confirm that there is an essential singularity located at r = 0. In addition, for large values of 
r, the Kretschmann and Ricci scalars yield 12vl 2 and 6A , respectively, which confirm that obtained 
solution is asymptotically AdS. In order to investigate the existence of horizon, one should consider 
vanishing metric function at the horizon. Considering / (r_|_) = 0, one can write 


to + Ar+ + 2 q 2 In 



= 0 . 


(47) 
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Fig. 5 Different scales of f(r) versus r for /} = 0.0001, l = 1 , A = —1 , q = 1 , and m = 0.1 (solid line), m = 1 
(dotted line), m = 3 (dashed line) and m = 8 (bold solid line). 




Fig. 6 Different scales of /(r) versus r for m = 3.1, Z = 1, A — —1, q = 1, and /3 = 0.001 (solid line), (3 = 0.003 
(dotted line), /3 = 0.006 (dashed line) and (3 = 0.009 (bold solid line). 


Obtaining analytical solutions of Eq. (H71) with respect to r + is not easy. So in order to investigate the 
horizons, we plot some graphs for the metric function (see Figs. [5] and [5]) • In addition, for studying 
the effects of the additional correction (/3—term), we plot Fig. [6] for metric function versus radial 
coordinate. 

As one can see, for small values of geometrical mass (Fig. [5]), the metric function has only one root 
which is located at small values of r. In general the effect of mass on number of horizon(s) could be 
divided into three regions that are denoted by two critical values for the mass parameter; m\ and m 2 
in which m\ < m 2 . For m = m\ and m = m 2 , two roots exist. In case of mi < m < m 2 , there are 
three roots for metric function. Otherwise, metric function has only one root. The location of outer 
(event) horizon is an increasing functions of geometrical mass. 

Next, due to our interest in nonlinear theory that was used in this paper, we plot the metric 
function diagrams for variation of the /3 (Fig. [H]) - It is evident from studying Fig. [(3 that the effect of 
the nonlinearity parameter is the same as geometrical mass and similar behavior is observed. 


















11 


4.1 Thermodynamics and conserved quantities 


Here, we are interested in studying thermodynamic properties of obtained black hole and calculating 
conserved quantities. First of all as one can see, the metric, that was employed for obtaining solution, 
contains a Killing vector field which is temporal (x M = <5 q ). Considering the concept of surface gravity, 
one finds the temperature of obtained black hole solution 




(V M X„) (Wx«) = 


r+A + r+q 2 + 2q 4 (3 
27rr^_ 


0(/3 2 ). 


(48) 


The next step is devoted to calculate the entropy of black holes. Since we are working in the Einstein 
gravity, it is allowed to employ the area law for calculating the entropy. By doing so, one can find the 
same relation such as that in Eq. (I18f) . It is notable that although entropy has the same form in both 
Maxwell theory and NLED, r+ in NLED is different from that in Maxwell one. Considering the flux of 
the electric field at infinity leads to similar relation for the electric charge as we obtained in Eq. (HU). 
It means that the nonlinearity does not affect the total electric charge of black hole. This behavior was 
also observed in other BI type theories of electrodynamics 1111- 

Next, we calculate the electric potential, at the horizon with respect to a reference with vanishing 
electric potential 

* = A^l™/e,ence - W\ r ^ r+ = ~Q (in (^) - + O (P) - (49) 

Following the same approach and applying the same counterterm action, we find that the nonlin¬ 
earity does not change the form of finite mass. Therefore, we obtain the finite mass as 

M = j, (50) 


where one can obtain m through Eq. S3- 

Now, we are in a position to study the validation of the obtained thermodynamic quantities, through 
the first law of thermodynamics. Using Eqs. (HU), HU, HU and considering M as a function of extensive 
parameters S and Q , we obtain 


M (S', Q) 


1 

2 


IAS 2 


+ 2Q 2 In 



2t t 2 Q 4 

~S 2 ~ 


P 


+ O (f3 2 


(51) 


where we considered /3, l and A as a fixed parameter. According to the first law of thermodynamics, 
we should check the following relation 


dM = TdS + <PdQ. 


(52) 


It is straightforward to calculate and y and check their results to be in agreement with 

what were previously obtained in Eqs. (1151) and (HU, respectively. Therefore, obtained conserved and 
thermodynamic quantities satisfy the first law of thermodynamics. 

According to the mentioned discussion for considering cosmological constant as a thermodynamical 
variable, one should generalize M (S, Q) to M (S, Q , A). Therefore Eqs. (1511) and (1551) will be modified 
as 


M (S,Q,A) = —- 


AS 2 „^ 2l (2 SV^A\ 2tt 2 Q 4 
— +2q 


+ o 03 2 ) . 


(53) 


and 


dM = TdS + <PdQ + OdA , 


d M\ 

~9sJq,a 


(p = 




where 


( 54 ) 
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Fig. 7 Different scales of 10 Cq and (bold lines) versus r+ for /3 = 0.001, A = —1, and q = 1 (solid line), 
q = 1.2 (dotted line) and q = 1.4 (dashed line). 


4.2 Thermal stability 


Using Eqs. (ED, ED and ED. one can find heat capacity of nonlinearly charged black hole solution 



In order to investigate the type one phase transition points, one should find the roots of the 
numerator of the heat capacity 

Ar A + q 2 r\ - 2q A f3 = 0. (56) 

Solving Eq. ED. one finds two real positive roots in which type one phase transition takes place 


r +lc Q =o = Q 


±X ~ 1 

2 A 


(57) 


where \ = + 8 A/3. Due to square root function and our interest in AdS spacetime, we have the 

restriction in which /? < —1/8 A. Next, as for the type two phase transitions, one can find roots of 
denominator of the heat capacity 

Ar% - q 2 r 2 + + 6q 4 /3 = 0, (58) 

which results into one real positive root 


r+\c Q - 


= q 


1 - A 
2 A ’ 


(59) 


where A = y/l — 24 Af3. Overall by what was mentioned, one can conclude that there are three phase 
transition points for these black holes: two of these phase transitions are type one and the other one is 
type two. We plot two diagrams for studying thermal stability in context of canonical ensemble (Figs. 
[T] and 0) • 

It is evident from studying the diagrams (Figs. [7] and [HD that in this case (NLED source), we are 
dealing with two types of phase transitions. One is related to the root(s) of the heat capacity and 
the other one is related to divergence points of the heat capacity. As one can see, there is a region 
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Fig. 8 Different scales of Cq and T (bold lines) versus r+ for q = 1, A = —1, and /3 = 0.001 (solid line), 
/3 = 0.003 (dotted line) and /3 = 0.006 (dashed line). 


(r+ < ri + where ri+ is the smaller root of the heat capacity) in which heat capacity is negative and 
system is unstable. In this region the temperature is positive. In r\ + < r. |_ < rd+ (where rd+ is 
divergence point), the heat capacity is positive but the temperature is negative. Therefore, the system 
is not in a real physical stable condition. Both the temperature and the heat capacity are negative 
in the interval r^ + < r_|_ < r%+ in which r 2 + is the larger root of heat capacity. Finally, the system 
acquires a physically thermal stable state in V 2 + < r+. In other words, in this region both heat capacity 
and temperature are positive and the system is in thermal stable state. It is notable that the divergence 
point of the heat capacity is located in the region in which temperature is negative. Therefore, it is 
not a physical phase transition. As one can see, in case of variation of electric charge (Fig. [TJ, the 
roots and the divergence points of the heat capacity are increasing functions of electric charge. For the 
nonlinearity parameter, the smaller root and divergence point are increasing functions of /3 (Fig. [Sj left) 
whereas the larger root is a decreasing function of (3 (Fig. [5] right). This shows that for increasing /3, 
the interval between larger root and divergence point (r^ < r+ < r 2 ) decreases and a compactification 
takes place. 

Next, we investigate determinant of Hessian matrix in the context of the grand canonical ensemble, 
in which finite mass of the black hole is a thermodynamical potential with two extensive parameters, 
S and Q. This assumption leads to following form of Hessian matrix 


H 


M _ 
S,Q — 


-AS 4 +Q 2 S 2 7r 2 -67r 4 Q 4 /3 2Q(4 tt 2 Q 2 /3-S 2 ) 

tt 2 S 4 S 5 

2Q(4ti- 2 Q 2 /3-S 2 ) S 2 ln(^)-127r 2 Q 2 /3 

S 3 S 2 


(60) 


which by using Eqs. (US and OH), we find 

I h m I OK - 9 2 4 + 6 <? 4 /3) , fr% \ V [(3 Arl + 5 q 2 ) p - r 2 + \ 

I h s,qI =-- ln {l^) + -^4-• (61) 

Using obtained determinant, one can plot its diagrams (Figs, [ill and [Toll in order to study thermal 
stability of the system. 

It is evident from studying these diagrams (Figs, [ill and flOl) . that there is a critical horizon radius, 
r +c , in which for r .(_ < r .|_ c the determinant of the Hessian matrix is positive which means the system 
is in thermal stable state. For the case of > r+ c , the determinant of the Hessian matrix is negative. 
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line), q = 1.2 (dotted line) and q = 1.4 (dashed line). 




Ih m I 

Fig. 10 Different scales of and 10T (bold lines) versus r+ for q M 1, l = 1, A = —1, and (3 = 0.001 

(solid line), /3 = 0.003 (dotted line) and /? = 0.006 (dashed line). 


r +c is an increasing function of the electric charge and nonlinearity parameter. It is worthwhile to 
mention that for a region of the r + , the temperature and determinant of the Hessian matrix are both 
positive which in that region the heat capacity is negative. As one can see in this case, thermal stability 
of the system in canonical and grand canonical ensembles does not have the same result. Therefore, 
we encounter a case of ensemble dependency. 

In order to remove ensemble dependency, we consider the assumption that was used in case of 
linearly charged BTZ black hole. One can find the mass of black hole as a function of three extensive 
parameters: electric charge, entropy and cosmological constant. Taking into account this assumption, 
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Fig. 11 Different scales of |H^q A | and T (bold lines) versus r+ for /3 = 0.001, A = —1, and q = 1 (solid 
line), q = 1.2 (dotted line) and q = 1.4 (dashed line). 


the Hessian matrix for nonlinear charged black hole will be as follow 


TT M 

n S,Q,A 


-AS i +Q 2 S 2 7v 2 -6n i Q i 0 2Q(4 k 2 Q 2 0-S 2 ) 

ir 2 S L S 3 


2Q(4ti- 2 Q 2 /3-S 

S 3 


') S 2 l„(^)-12. 2 Q 2 /3 

S 2 


Q 

A 


Q 
' A 


Q 
2 A' 2 J 


where its determinant is 


IttM I 

| Jrl S',Q,yl| 


A In (—+ 2 B 
87T 2 A 2 r 4 


in which 


(62) 


(63) 


A = 2 A 2 r% - q 2 (,fr\ - Ar% - 6g 4 /3), 

B = 16 q 6 (3 + (22 Afi - 3) q A r\ + 3 (4H/3 - 1) Aq 2 r\. 

Now, we plot the related diagrams of Eq. (TlITIfl and compare them with the results of the heat 
capacity (see Figs. [T][9]). 

As one can see, considering this modification changes the number of the roots of the determinant 
of the Hessian matrix. In other words, contrary to previous case (Figs. [9l and fTOT) in which the deter¬ 
minant of the Hessian matrix had only one root, in this case it has two roots (Figs. fiTmi) . We name 
these two roots r\+ and r 2 + in which ri+ < r 2 +. For n+ < r+ < r 2 +, determinant of the Hessian 
matrix is negative, otherwise (r + < r\+ and r+ > r 2 +) it is positive. Interestingly, by considering cos¬ 
mological constant as a thermodynamical variable, the larger root of the heat capacity, temperature 
and determinant of the Hessian matrix coincide. In other words, for r + > r 2 + in both canonical and 
grand canonical ensembles, the system is thermally stable and temperature is positive. Therefore, the 
ensemble dependency is removed by this consideration. 

Considering the cosmological constant also enable us to consider the variation of it in plotting 
diagrams. As one can see in this case, the smaller root of the determinant of the Hessian matrix and 
temperature are decreasing functions of the cosmological constant (Fig. [13l left). On the other hand, 
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Fig. 12 Different scales of |H^q A | and T (bold lines) versus r+ for q = 1, A = — 1, and p = 0.001 (solid 
line), p = 0.003 (dotted line) and p = 0.006 (dashed line). 




Fig. 13 Different scales of and T (bold lines) versus r+ for q = 1, P = 0.005, and A = —1 (solid 

line), A = —2 (dotted line) and A = —3 (dashed line). 


the larger root of the determinant of the Hessian matrix and temperature of the system are increasing 
functions of A. 

Here, we discuss the effects of considering the nonlinearity parameter, p, as a thermodynamical 
variable. In other words, taking into account the Hessian matrix as a function of /?, one can regard 
H = H(S, Q, P) or H = H(S, Q, A, p). Considering (5, Q, P) and ( S , Q, A, /?) as extensive parameters, 
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Fig. 14 Different scales of 10 k x |Hg[q j/3 | and T (bold lines) versus ?'+ for /3 = 0.001, l = 1, A = —1, and 
q = 1 (solid line), q = 1.2 (dotted line) and q = 1.4 (dashed line). ”In order to plot clear figures, we set k = —1 
and k = — 5 for left and right figures, respectively” 


one finds their related hessian matrices are, respectively, 


IT M 

n S,0,? 


-4S 4 +Q 2 S 2 7r 2 -67r 4 Q 4 /3 2Q(4 ti 2 Q 2 /3-S 2 ) 


jt 2 S 4 

2Q(4ti 2 Q 2 /3-S 2 ) 

S 3 


TS 3 " 


S 3 


In 


(#) 


127r 2 Q 2 /3 


S 2 


4Q 3 tt 2 

S 3- 


2Q 4 tt 2 

S 3 

4Q 3 tt 2 

S 2 

0 


(64) 


and 


H 


M 

S,Q,P,A 


-ylS 4 +Q 2 S 2 7i 2 -67r 4 Q 4 f 

2Q(4tt 2 Q 2 /3-S 2 ) 

-s 3 - 


2Q 4 tt 2 


S 3 
S 


Straightforward calculations show that 


I tt.M I 


AAq 6 

TT 2 r\ 


In 


2Q(4tt 2 Q 2 /3-S 2 ) 

S 3 


7 —7T 2 \ _ 12tt 2 Q 2 /3 


) 


S 3 


4Q 3 7r 2 

S 2 


Q 

A 


In 


(*) 


+ 4 


7r 2 r^ 


2Q 4 7T 2 

“s 3- 

4Q 3 tt 2 

s 3- 

0 

0 


4g 10 /3 

7r 2 r® 


Q 
' A 


<3 

2 A 2 


(65) 


( 66 ) 


IttM I 
|- tl S,Q,A,/3| 


3 q 8 


2„4 


7r^r 


2yl7r 2 r^ 


^ g 10 [In (—r+Tl) + 6] ^ g 12 ,0 


87l7r 2 r^_ 


2An 2 r s f ’ 


(67) 


Now, we plot Figs, fill and fTKl to discuss ensemble dependency. Regarding Figs. H4l and fl5l we find 
that considering /3 as a thermodynamical variable (with or without regarding A as a thermodynamical 
parameter) leads to ensemble dependency. In other words, we conclude that /? is not a thermodynamical 
variable. 
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Fig. 15 Different scales of 10 fc x |Hg f Q A and T (bold lines) versus r+ for /3 = 0.001, A = —1, and q = 1 
(solid line), q = 1.2 (dotted line) and q = 1.4 (dashed line). ”In order to plot clear figures, we set k = —1 and 
k = — 5 for left and right figures, respectively” 


5 Summary and Conclusion 

In this paper, we have investigated thermodynamic properties of the charged BTZ black hole and 
studied thermal stability of it in the context of both canonical and grand canonical ensembles. It was 
observed that this black hole has two types of the phase transition; one was related to the changing 
in the signature of the heat capacity (vanishing heat capacity). There was also another type of the 
phase transition which was related to the divergency of the heat capacity. It was shown analytically 
that there exists a real positive root and also, a divergence point for the heat capacity. 

One of the most interesting results of this paper was ensemble dependency of thermal stability. It 
was seen that considering the total mass of black hole as a function of two extensive parameters (S & 
Q ) will lead to a 2 x 2 Hessian matrix. The determinant of this matrix was positive only for a region 
of small values of horizon radius. In other words, the conditions for thermal stability of the black hole 
for these two ensembles were different. Therefore, there was a case of ensemble dependency. In order to 
solve this ensemble dependency, we considered cosmological constant as an extensive thermodynamical 
parameter. This lead to the Hessian matrix being 3x3. This consideration resulted into removing 
ensemble dependency and both ensembles yield the same result. In other words, in both pictures, 
thermally stable physical black hole only observed for large r + . It is notable that the constant l in the 
gauge potential and logarithmic part of metric function inserted by hand for the reason of obtaining 
dimensionless argument for the logarithmic function. Calculations showed that this length parameter 
should be related to cosmological constant (l = l/y/—A) to remove ensemble dependency. Another 
important property is that our study support the idea that the cosmological constant is not a fixed 
parameter, but an extensive thermodynamical variable. In other words, the variation of this parameter, 
should be taken into account in studying thermodynamical behavior of the system. 

In the next part of the paper, we generalized Maxwell theory to the case of NLED. The Lagrangian 
of the mentioned nonlinear model was a quadratic Maxwell invariant term in addition to the Maxwell 
Lagrangian. It was seen that the structure of the metric function, number of the horizons and the 
behavior of the metric function were modified. In context of thermodynamical quantities, no effect on 
the forms of total mass, electric charge and entropy was observed whereas the temperature and electric 
potential were modified. 

In case of the phase transition in AdS spacetime it was seen that: BTZ chargeless black holes 
have no phase transitions. Linearly charged ones only have a type one phase transition. Whereas, 
nonlinearly charged black holes enjoy two types of phase transitions. On the other hand, being dS and 
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AdS spacetime highly modified the number and the type of phase transitions. This emphasizes the 
contributions and effects of background spacetime on thermodynamical structure and behavior of the 
system. 

In studying the stability, similar to linearly charged BTZ black holes, ensemble dependency was 
observed for nonlinearly charged BTZ solutions. In order to remove this ensemble dependency, we 
employed the method that was used in case of linearly charged BTZ black holes. It was seen that this 
method successfully removed ensemble dependency. This shows that our consideration of cosmological 
constant as thermodynamical variable is an acceptable one. In other words, cosmological constant 
indeed is a thermodynamical variable, not a fixed parameter and its variation should be taken into 
account in the first law of thermodynamics 

dM = TdS + <PdQ + OdA , 


where 


0 = 



We also checked the effects of considering the nonlinearity parameter, /3, as a thermodynamical variable. 
We found that regarding /3 as a thermodynamic (extensive) parameter leads to existence of ensemble 
dependency. Therefore in order to remove ensemble dependency, we should consider ( S , Q , A) as a set 
of extensive parameters. 
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